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'^The  effect  of  randomly  missed  range  data  on  target  tracking  performance 
has  been  investigated.  The  specific  problem  addressed  concerns  SEAFIRE.  Most 
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other  purposes. 
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r 


DD  1473  EDITION  OF  I NOV  69  IS  OBSOLETE 

S/N  0|02*0l4-660l 


UNCLASSIFIED 


security  CLASSIFICATION  OF  THIS  RAOt  fWhan  Dafa  Bniarad) 


7ft 


-I-V-IJMITY  CL  ASSiriC*  rcON  of  This  PAOerWA^n  Pmlm  en(««rf; 


20.  (Continued) 

the  target.  Both  of  these  models  provide  a quite  general  description  of  the 
stochastic  processes  being  analyzed.  The  effect  of  uniformly  (non-random) 
missed  range  data  was  also  explored  in  order  to  provide  a reference  which  is 
somewhat  easier  to  analyze  and  understand.  For  the  sake  of  completeness, 
angle  tracking  performance  was  briefly  examined  for  the  case  of  uninterrupted 
angle  measurements. 

Results  presented  in  this  technical  report  are  designed  to  aid  in 
determining  performance  to  be  expected  from  a target  tracking  algorithm 
subjected  to  randomly  interrupted  range  data  and  to  aid  in  specifying  sensor 
accuracy  requirements  necessary  to  obtain  the  desired  tracking  accuracy. 
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I. 


INTRODUCTION 


The  primary  purpose  of  a tracking  filter  in  any  fire  control  system 
is  to  process  sensor  information  to  obtain  target  state  estimates 
(position,  velocity,  etc.)  suitable  for  accurately  predicting  future 
target  position.  Designing  a tracking  filter/predictor  algorithm  is 
fairly  straight-forward  once  a target  model  and  sensor  error  model  have 
been  established.  The  above  is  true  even  for  the  case  of  randomly 
interrupted  range  data.  Generally,  the  tracking  filter  is  designed  to 
simply  "coast"  during  the  data  interruption.  Error  analysis  for  such  a 
system,  however,  is  quite  difficult  and  is  the  major  topic  addressed  in 
this  technical  report. 

For  the  case  of  uninterrupted  data  (assuming  the  usual  linear  target 
and  measurement  model  and  linear  estimation  algorithm) , filter/predictor 
error  covariance  (a  matrix  describing  the  random  nature  of  filter/ 
predictor  estimates)  is  deterministic  and  can  be  easily  calculated.  For 
the  case  of  randomly  interrupted  data,  however,  the  error  covariance  is 
not  determinisitc  but  a random  process.  It  is  precisely  this  property 
that  makes  error  analysis  of  the  interrupted  data  system  difficult. 

The  analysis  presented  in  this  report  computes  the  first  moment 
(expected  value)  of  the  filter/predictor  error  covariance  using  an  exact 
and  an  approximate  technique.  The  approximate  technique  is  necessary 
because  of  the  tremendous  computational  burden  associated  with  the  exact 
technique. 

Most  of  the  analysis  presented  is  directed  toward  range  filtering 
with  randomly  interrupted  range  data.  A brief  analysis,  however,  is  pre- 
sented for  angle  filtering  with  uninterrupted  angle  data. 
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11. 


TAKuKT  MOUtL 


Tht*  trtrijt't  nKKiol  .ii!sum»Hl  in  this  stiuiy  is  >■»  ixirrolatinl  .icoolt>r.tt  ion 
model  developed  for  tdiijot  tr.ickiiui  .analysis  by  Sinqer  in  l‘>70.  This 
model  was  chosen  K'cavjse  of  its  simplicity,  flexibility  and  previous  suc- 
cessful filter  design  applications  such  as  tlu'  d.ir  filter  (refer«'nce  2). 
An  uncouplevl  spherical  cvxirdiiiate  system  was  chosen  also  iH'cause  of 
simplicity  and  previous  successful  application  (reference  1).  As  an 
example  of  this  technique,  tlie  range  direction  mi'deling  assumptions  are 
illustrated  by  the  following  equation. 


(1) 


where  r is  t.irget  r.»nqe,  u'(t)  is  white  ikaussi.ui  process,  i is  .»  seri.nl 

"time  const.nnt"  and  superscript  .lenotes  time  dt'rivative. 

St.ntistic.il  properties  of  the  alv>ve  process  .ire  .-onvenit'nt  ly 
described  by  the  autocorrelation  function 


where  K(-)  doiiotes  statistic.nl  I'xpected  v.ilue  operation.  The  p.nr.imeter 
of  equation  (2)  specifies  the  target  "nvineuver  level"  (references  l,.'l. 

Modeling  for  the  elevation  .uui  lH*aring  directions  is  identical  to 
range  direction  mixleling  except  for  slight  differences  which  .nre  described 
in  reference  1. 


III.  MODKLINi.;  'HIK  MISSKD  DATA  PROCESS 
A first  order  Miirkov  process  was  chosen  to  model  the  randomness  of 
the  missing  data  process.  Such  a process  is  reasonably  versatile  in  that 
it  allows  the  probability  of  receiving  a measurement  and  correlation  time 
to  be  chosen  independently.  Also,  by  using  the  Chapman-Kolmogorof f equa- 
tion as  described  in  reference  4,  propagation  of  measurement  probabilities 
can  be  conveniently  described  by 


p(m^) 


P (m^) 


where  p(m^)  denotes  probability  of  receiving  a measurement  at  time  t^  and 
p(m^)  denotes  the  probability  of  not  receiving  a measurement  at  time  t^. 
The  2x2  transition  m.itrix  is  given  by 


p("V^iK) 


where  the  bar  "|"  indicates  a conditional  probability,  i.e.,  P 
denotes  probability  of  not  receiving  a measurement  at  time  given  a 

measurement  was  received  at  time  t^.  Appendix  A describes  how  the  elements 
of  X are  chosen  to  match  the  properties  of  a physical  system. 
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IV.  TRACKING  (I’lLTER/PKEDICTOR)  ALGORITHM 


This  study  assumes  a Kalman  tracking  filter  designed  about  the  un- 
coupled spherical  target  modeling  scheme  discussed  earlier.  Figure  1 
illustrates  the  structure  of  such  a filtering  technique  via  block  diagram. 
The  indicated  coupling  between  filters  is  included  to  allow  constant  maneu- 
verability modeling  as  discussed  in  reference  1.  The  coupling  is  quite 
"loose",  however,  as  demonstrated  by  Quigley  (reference  3)  where  a suc- 

! cessful  tracking  algorithm  was  designed  by  completely  ignoring  the  coupling. 

\ 

t 

[ Also,  it  is  assumed  the  prediction  function  is  performed  via  the 

I 

^ simple  state  transition  method.  This  technique  (in  spherical  coordinates) 

was  demonstrated  to  be  substantially  superior  (for  point  defense  applica- 
[ tions  (reference  3))  to  similar  techniques  in  Cartesian  coordinates. 

f 

The  analysis  presented  in  this  note  pertains  specifically  to  tlie 

i 

filter/predictor  assumptions  described  above.  The  results,  however,  are 
quite  indicative  of  performance  to  be  expected  from  other  well  designed 
tracking  schemes  addressing  the  same  problem. 
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measured  bearing 
at  time  tj^ 


Bearing  Coordinate 
Kalman  Filter 


measured  range 
at  time  tj^ 


Range  Coordinate 
Kalman  Filter 
("coast"  during 
missing  data) 


measured  elevation 
at  time  tj^ 


Elevation  Coordinate 
Kalman  Filter 


Figure  1.  Kalman  Tracking  Algoritlim  Using  Uncoupled 
Spherical  Model  (Block  Diagram) 
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V.  EXACT  SOLUTION  KOK  ERkOK  COVARIANCE  WITH  RANDOMLY  MISSED  DATA 

The  error  covariance  matrix  of  a Kalman  filter  algorithm  is  propagated 
(in  time)  by  the  following  equations: 

• Vl  “kU  <“k*l  "k*l  "kkl  * 

■’k.i  ■ » - \*i  "kti'  "k« 

(1  = identity  matrix) 

The  filter  error  covariance,  P,  is  defined  by 

Pu  = k[  (^  - ^)  (^  - )'^] 

k *■  — k — k -k 

where 

= system  state  vector  (position,  velocity,  etc.)  at  time  ij 

' ■) 

X,  = filter  estimate  of  X at  time  t,  . |J 

— k — k 

r i 

EL ’J  denotes  expected  value  operation.  5 

Superscript  "T"  denotes  transpose.  |j 

• 1 

The  matrix  quantities  R and  Q are  the  system  state  transition  (j 


matrix,  measurement  error  covariance  matrix  and  process  disturbance  co- 
variance  matrix  respectively. 

Definitions  and  physical  significance  of  cji , R and  Q are  described 
in  references  5 and  6 while  Singer  (reference  1)  describes  the  exact 
formulation  of  these  parameters  for  the  target  model  being  considered. 
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The  matrix  H of  equations  (6)  is  the  observation  matrix  which 
indicates  how  measurements  are  related  to  the  system  state  vector.  As 

( r ) 

an  example,  consider  that  noisy  measurements  of  range  (Z  ) are  made 
and  the  range  coordinate  state  vector  consists  of  range  (r) , range  rate 
(r)  and  range  acceleration  (r) . The  measurement  equation  becomes 


k+1  k+1  -k+1 


= [l  I 0 I o] 


(7) 


where  represents  the  error  in  measured  range  at  time  . Notice 

that  if  the  measurement  at  time  is  missed,  becomes  all  zero  and 

equations  (6)  reduce  to 


"k.l  = " ^k- 


For  the  case  of  uninterrupted  data  (or  non-rai.^omly  interrupted  data) 
the  filter  error  covariance  is  a deterministic  quantity  fully  described  by 
equations  (6) . For  the  case  of  randomly  interrupted  data,  however,  the 
situation  is  considerably  more  complex  since  error  covariance  becomes  a 
random  process.  An  exact  solution  for  the  expected  value  of  filter  error 
covariance  is  easily  formulated  as  the  weighted  average  of  all  possible 
filter  error  covariances.  Such  a formulation  is  given  by 


^ Prob(case.) 

all  cases 


(9) 
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where  case^  denotes  the  possible  measurement  sequence  and  Prob(case^) 
denotes  the  probability  of  case^  occurring.  To  clarify  the  above,  consider 
the  simple  excimple  where  is  known  cind  the  expected  value  of  P^  is 
desired.  The  following  two  measurement  sequences  (cases)  are  possible. 


Case^^  measurement  does  not  occur  at  time  t^. 

Case2  measurement  does  occur  at  time  t^. 

From  equations  (6),  Pj^(case^)  for  i = 1,2  can  be  formulated  to  yield 


T 

P^(case^)  = <(>Pq4>  + 

P^Ccase^)  = (I  - 
T 

where  M,  = <{ip  <)>  +0 

1 o o 

Ki  = 

Assuming  the  probability  of  missing  a measurement  at  time  t^  is  q^  (which 
implies  the  probability  of  receiving  a measurement  at  time  t^  is  1 “ <3^^ 
results  in  the  following,  exact  expression  for  the  expected  value  of  P^^. 


i 


1 


eCp^]  = (jl-q^)  (I-K^H^)M^]  + q^[<tP^<Ji'^  + ] 


where  and  are  defined  earlier.  Prom  this  simple  example,  it  is  easy 
to  see  the  difficulty  in  implementing  equation  (9)  for  k fairly  large. 
Specifically,  the  complexity  of  equation  (9)  increases  exponentially  with 
increasing  time  steps.  For  the  general  case. 


= r 


. ^k+l 
i=2 


k+l-*  (case^)  Prob(case^). 
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Prediction 


The  prediction  error  covariance  equations  are  linear  and  therefore 
allow  for  a simple  solution  for  the  expected  value.  Specifically,  the 
expected  value  of  prediction  error  covariance  is  given  by 

E[M(tp)]  = <C(tp)E[P^]<t'’'(tp)  + Q(tp) 

where 

M(tp)  =E[(X(t^+tp)  -X(VV^‘^‘VV 

t = prediction  time 

P 

t = time  of  last  measurement 

k 

Restrictions 

In  practice,  it  has  been  found  that  equation  (9)  has  practical  use 
for  k * 10  or  less.  It  is  precisely  this  computational  restriction  that 
motivated  development  of  the  recursive  approximate  solution  presented  in 


the  next  section. 


VI.  RKCURSIVli  APPROXIMATE  SOLUTION  EOR  ERROR 
COVARIANCE  WITH  RANDOMLY  MISSED  DATA 
The  covariance  propagation  equations  (equations  (6))  can  bt’  comluned 

to  yield  one  nonlinear  equation  for  the  case  of  measurement  occurring  at 

time  • This  equation  is  shown  below  in  an  explicit  and  an  abbreviated 

form. 


(11) 


Vl  = 'NL<"k) 


(12) 


where  f (•)  represents  the  nonlinear  function  described  in  e^niation  (11). 

NL 

From  equation  (8)  (case  of  no  measurement  at  time  wt'  liave  tlie  linear 

equation 

'kki  ■ 

’’k.i  ■ 

where  f.  (•)  represents  the  linear  function  described  in  equation  (13).  Using 

Li 

the  notation  jest  defined,  the  expected  value  of  can  be  expressed  by 

E[P,^^l]  - e[  f^^(P^)]  prob(m^^^)  + E[fj^(P^)  ] Prob(m^^^)  (15) 
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where  prob(mj^^j^)  denotes  the  probability  of  receiving  a measurement  at  time 
tj^^^  and  prob(m^^j^)  denotes  the  probability  of  not  receiving  a measurement 
at  time  To  generate  the  recursive  approximate  solution,  the  following 

approximation  is  needed. 

This  approximation  allows  equation  (15)  to  be  rewritten  yielding  the  very 
simple  recursive  solution 


(16) 


Prediction 

Expected  value  of  the  prediction  error  covariance  is  obtained  exactly 
as  for  the  exact  solution  of  the  previous  section. 


Restrictions 

This  approximate  solution  effectively  ignores  2nd  and  higher  order 
statistical  moments  of  the  matrix  P.  Appendix  B presents  experiments 
conducted  to  determine  accuracy  of  the  recursive  approximate  solution 
(equation  (16)).  From  these  experiments,  it  was  concluded  the  approxi- 
mate solution  is  quite  accurate  enough  for  the  intended  purposes. 
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VII.  RESULTS 


This  section  primarily  seeks  to  describe  the  effect  of  missing  range 
data  on  tracking  filter /predictor  performance.  A very  brief  analysis  of 
angle  tracking  is  also  presented.  The  cases  of  randomly  missed  and  uni- 
formly missed  range  data  are  explored.  Experiments  were  conducted  in- 
volving surface  and  air  target  models.  Both  models  were  Singer  (reference  1) 
correlated  acceleration  models  with  the  following  parameters: 


Air 


0 =0.5  yds/sec 
m 

1 = 20  sec 
m 


Surface 


2 

0 =0.1  yds/sec 

m 

T =20  sec. 
m 


These  parameters  were  chosen  based  on  experience  gained  during  the  Gunnery 
Improvement  Program  (GIP)  (reference  2) . For  both  models,  the  uninter- 
rupted data  rate  was  assumed  16  Hz.  The  remainder  of  this  section  presents 
graphical  summaries  of  the  experiments  conducted. 

Randomly  Interrupted  Range  Data 

Results  of  these  experiments  were  obtained  using  steady  state  solu- 
tions obtained  via  the  recursive  approximate  technique  described  in  the 
previous  section  (equation  (16))  and  in  Appendix  B.  The  results  are 
presented  by  standard  deviation  prediction  error  contours.  These  contours 
present  the  data  in  a very  compact  form  and  therefore  must  be  interpreted 
very  carefully.  Each  contour  represents  a locus  of  necessary  conditions 
for  a particular  performance  criterion  (standard  deviation  prediction 
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error  in  position) . A contour  may  be  interpreted  as  a plot  of  required 


i 


probability  of  receiving  a measurement  as  a function  of  sensor  error 
standard  deviation.  Also,  a contour  may  be  interpreted  as  a plot  of 
required  standard  deviation  sensor  error  as  a function  of  probability 
of  receiving  a measurement.  Specific  procedures  for  interpreting  these 
contours  are  best  illustrated  by  example. 

Example  1;  Refer  to  Figure  2 which  addresses  a prediction  time  of 
25  sec  for  surface  targets.  Assume  the  sensor  of  interest  provides  range 
data  with  errors  of  27.5  yds  standard  deviation.  Locating  27.5  yds  on  the 
horizontal  axis  and  projecting  vertically  results  in  intersecting  the  40 
yd  contour  at  prob (measurement)  = 0.7.  This  result  indicates  that  the 
given  sensor  must  provide  probability  of  measurement  ^0.7  to  achieve  an 
expected  standard  deviation  prediction  error  of  40  yds.  Also,  notice  that 
the  30  yd  contour  is  not  intersected  for  the  sensor  of  interest.  This 
result  indicates  that  such  performance  is  not  possible  for  the  given  sensor 
even  if  all  measurements  are  received. 

Example  2:  Refer  to  Figure  3 which  addresses  prediction  time  of 
10  sec  for  air  targets.  Assume  the  sensor  of  interest  provides  a proba- 
bility of  measurement  = 0.0,  Locating  0.9  on  the  vertical  axis  and  pro- 
jecting horizontally  results  in  intersecting  the  20,  30  and  40  yd  contours 
at  3.5,  17.5  and  42  yds  standard  deviation  sensor  error.  These  results 
indicate  that  achieving  standard  deviation  prediction  errors  of  20,  30  and 
40  yds  requires  standard  deviation  sensor  errors  at  least  as  good  as  3.5, 
17.5  and  42  yds  respectively. 
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, , , 1 T 1 — 

5 10  15  20  25  30  35 

STANDARD  DEVIATION  OF  MEASUREMENT  ERROR  (YDS) 


Kiqure  2.  Standard  Deviation  Prediction  Krror  Contours 
for  Surface  Target  (25  sec  prediction  time) 
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Unifoimly  Intorrupteti  Kjiujo  tota 

Kosults  of  tlu'so  oxpor imt'iits  wero  obtaiiuni  usimi  stoady  state  soU»- 
t ions  obtained  from  the  standard  K.tlman  filter  covariance  equations 
(equations  (n)  ) . The  time  step  lx.‘tween  samples  was  varied  from  1 'It-  sec 
to  1 sec  to  represent  uniformly  missinq  tu'  data  to  uniformly  missinu  IS 
of  lo  samples.  Hotli  surface  and  air  taruets  were  t.\>ns idereii  for  variovis 
prediction  timt's.  Results  are  displayed  in  I’iqures  4 (surface)  and  *'  (air) 
wl\ere  standard  deviations  of  predict iot>  error  are  sliown  as  a function  of 
time  step  interval  for  various  prediction  times  and  sensor  accuracies. 
Aqain,  t)ie  procedure  for  interpret  inq  file  qr.>pl>s  will  Iv  illustrated  by 
example. 

K.x.tmple:  Refer  to  I'iqure  S wlucli  addresses  r.inqe  traclvinq  of  .nr 

tarqets  with  uniformly  missed  ranqe  data.  Assvune  tli.vt  d.»ta  is  unifojtnly 
received  every  .1 '4  sec.  Kurtlier,  assiune  .i  prediction  time  of  10  sec  and 
sensor  st.tndard  di'viation  error  of  S yds.  U'c.it  inq  -) '4  sec  on  the  .ibscissa 
and  projectinvi  vertic.illy  to  intersect  the  10  sec,  f>  yd  st.ind.ird  di'viation 
error  curve  reve.ils  st.ind.ird  deviation  prediction  error  is  .10  yvis. 

Sensitivity  Analysis 

A brief  sensitivity  analysis  was  performed  to  determine  t)ie  effect 

of  r and  c of  Sinqer's  tarqet  model  (reference  1)  v'n  ste.ulv  state  pre- 
m m ■ 

diction  error  vari.ince.  The  main  objective  of  tliis  an.ilysis  was  to  ex.imine 

sensitivity  witli  respect  to  since  this  par.imeter  is  the  nx^st  .lifficult 

to  specify.  Kor  convenience,  the  uniform  missed  d.tta  solution  w.is  used. 

Results  vire  presented  in  Kiqures  e throuqh  where  V'erformance  qraphs,  sue)! 

as  just  discussed,  .ire  presented  for  various  values  of  i and  . Prom 

m m 

these  results  it  was  concludixl  tli.it  sensitivity  to  r 

In 


ij 


i 

t 


is  relatively  mild. 


STANDARD  DEVIATION  PREDICTION  ERROR  (YDS) 


STANDARD  DEVIATION  PREDICTION  ERROR  (YDS) 


STANDARD  DEVIATION  PREDICTION  ERROR  lYDS) 


STANDARD  DEVIATION  PREDICTION  ERROR  (YDS) 


Bearing  Direction  Analysis  (uninterrupted  data) 


For  the  sake  of  completeness,  a brief  analysis  was  conducted  to  j 

i 

determine  bearing  direction  prediction  accuracy.  The  analysis  was  performed 

for  surface  mode  only  since  an  air  target  analysis  would  require  assuming  a 1 

I 

specific  target  trajectory.  Results  are  shown  in  Figure  10  where  standard 
deviation  of  prediction  error  is  plotted  against  prediction  time  for  various 
sensor  statistics.  The  abscissa  also  displays  an  intercept  range  scale 

which  was  obtained  from  reference  7.  An  objective  of  this  analysis  was  to  ! 

determine  a practical  lower  limit  for  sensor  angular  accuracy,  i.e.,  find  | 

an  accuracy  such  that  better  accuracy  does  not  improve  predictor  performance. 

From  Figure  10  it  can  be  seen  that  such  an  accuracy  is  quite  stringent  (at  ^ 

least  below  1/8  millirandian) . Since  this  accuracy  is  probably  lower  than  | 

1 

practical  to  achieve,  no  extended  effort  was  made  to  find  the  actual  lower 
limit. 
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standard  deviation  prediction  error  (MILLIRADIANS) 
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1 \ »Pr«diction  Tim*  U*c) 

2S  30 

-I » f ^ Rang*  (Kyd*) 

12  13  14 


Fiqure  10.  Bearing  Direction  Brodiction  Accuracy  (surface  target) 


Vlll.  CONCUlUINO  R1>1ARKS 


Several  techniques  have  been  developed  for  analyzing  the  effects 
of  missing  range  data  on  tracking  filter/predictor  performance.  These 
techniques  were  exercised  via  the  CDC  6700  computer  to  provide  performance 
estimates  for  the  particular  conditions  investigated.  Computer  programs 
for  estimating  filter/predictor  perfornutnce  for  other  conditions  of  interest 
are  available  on  request. 
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APPENDIX  A 


RELATING  TRANSITION  PROBABILITIES 
TO  A PHYSICAL  SYSTEM 


The  purpose  of  this  appendix  is  to  relate  the  transition  probabilities 
of  equation  (5)  to  parameters  describing  a physical  system.  The  physical 
system  selected  to  describe  the  missing  data  process  is  shown  in  Figure  A-1. 

Basically,  the  selected  system  consists  of  a zero  mean  unit  variance  1st  i 

order  Gauss  Markov  process  driving  a threshold  detector  which  determines 
if  a measurement  is  available  at  a particular  time  step.  The  1st  order 
Gauss  Markov  process  is  conveniently  described  by  its  autocorrelation 
function 

From  the  autocorrelation  function,  ^ can  be  specified  in  terms  of  an  auto- 
correlation time  constant,  r,  by 

/ Ml 

y-  = c“  r 

where  At  is  the  sample  time  interval  ~ 

The  specific  objective  of  this  appendix  is  to  relate  the  matrix  T of 
equation  (5)  to  the  physical  system  jiarameters  i (autocorrelation  time 
constant)  and  T (threshold  value) . 


A-1 


Solving  for  T 


The  upper  left  element  of  T (equation  (5))  is  defined  by 


I'll  ^ prob{nv^^Jm^}=  prob{x^^^>T  | x^>T) 


where  prob{a|b}  denotes  the  probability  that  event  "a"  occurs  given 
event  "b"  occurred.  The  expression  (A-1)  can  be  expanded  via  Bayes'  rule 
to  yield 


prob{  Xj^^l>T,  Xj^>T} 
11  prob{x^>T] 


where  prob{a,b}  denotes  the  probability  of  event  "a"  and  "b"  both  occurring. 
The  denominator  term  of  equation  {A-2)  is  simply  the  unconditional  proba- 
bility of  receiving  a measurement  and  is  given  by 


prob{Xj^>T}  = / {0,l)dXj^ 


X,  =T  k 
k 


where 


Applying  Bayes'  rule  again,  the  numerator  term  of  equation  (A-2)  eventually 
becomes 

V”  \+i"~ 

prob{Xk^l>T,  x^>T}  - / N ^ 

\+l=T 


2 2 

where  o = 1 - ^ . Combining  equations  (A-4) , (A-3)  and  (A-2)  yeilds  the 
final  expression 


/ N (0,1)  / 


N (fx^,  a ) d d 

X,  , ' k (1)  X,  , X, 

Xj^=T  k *k+l~^  ^ 


11 


X 

^ N^(0.1)d^ 

x^=T 


The  parameters  ^ and  o are  related  to  the  physical  parameters  by 


2 AT 

T 


o = 1 - e 
(0 

AT 

f = e T 


and  At  is  sample  interval  time  step. 


^22'  ^12'  ’’^21 


It  is  easily  shown  that  the  remaining  elements  of  T are  given  by 


T = 1 - T 
21  11 


12 


(1-T^j^)  prob(x>T) 
l-prob(x>T) 


T = 1 - T 
22  12 


prob(x>T)  = / N (0,l)dx 
x=T 


A-4 


Thus  the  proUibility  transition  matrix  I'  has  U't'n  oxprosse.i  in  torms  of 
the  physical  systt'ra  threshold  value  (T)  and  autocorrelation  time  constant 
(t). 
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APPENDIX  B 

ACCURACY  OF 

RECURSIVE  APPROXIMATE  SOLUTION 


The  purpose  of  this  appendix  is  to  provide  an  accuracy  evaluation  of 


the  recursive  approximate  solution  described  by  equation  (16) . it  is 
easily  shown  that  the  recursive  approximate  solution  (RAS)  is  exact  for 
the  cases  of  unity  or  zero  probability  of  receiving  measurements.  For 
other  cases,  however,  accuracy  evaluation  is  much  more  difficult.  The 
evaluation  technique  described  in  this  appendix  consists  of  comparing 
the  RAS  with  the  exact  solution  of  equation  (9)  for  early  transient  and 
steady  state  responses.  The  comparisons  presented  here  are  based  on  the 
following  modeling  assumptions. 

TARGET  MODEL 

Singer  correlated  acceleration  model  (see  reference  1)  with 

2 

= 0.1  yds/sec  (maneuver  level) 

= 20  sec  (maneuver  correlation  time  constant) 

The  measurement  sequence  is  at  16  Hz  with  standard  deviation  error  of 
10  yds. 

Experiment  #1  (Early  Transient  Response)  i; 

This  experiment  serves  to  evaluate  the  early  transient  response  of 
the  recursive  approximate  solution.  Initial  state  error  covariance  is 
given  by 

f 

j 


B-1 


1 

0 0 

2 

(30  yds/sec)  0 , 

0 (0.1  yds/sec*”)^ 

The  miss  data  process  autocorrelation  time  constant  (see  Appendix  A) 
is  10  seconds.  The  evaluation  presented  here  consists  of  a comparison 
between  the  recursive  approximate  solution  (equation  (16))  and  the  exact 
solution  equation  (9).  Expected  standard  deviation  of  filter  error  in 
position  (o^)  and  velocity  (o^.)  for  the  two  solutions  is  plotted  in 
Figures  B-1  and  B-2.  Expected  standard  deviation  of  acceleration  error 
is  virtually  identical  for  both  solutions  and  was  therefore  not  plotted. 

Figure  B-1  results  from  prob (measurement)  =0.9  while  Figure  B-2  results 
from  prob (measurement)  = 0.5.  Note  that  in  both  cases  the  difference 
between  the  two  solutions  is  small. 

Experiment  #2  (steady  state  response) 

Evaluating  accuracy  of  the  steady  state  solution  is  a bit  more  compli- 
cated. This  complication  arises  because  the  following  principle  cannot 
be  easily  deduced  from  the  given  exact  formation  (equation  (9))  for 
finite  jc. 

\ 

where  is  the  filter  error  covariance  matrix  at  time  step  )c  and  £[•] 
denotes  expected  value. 
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Thus  recognizing  a steady  state  solution  (for  finite  K)  using  the 
formUjlation  of  equation  (9)  is  impossible. 

The  technique  used  here  for  testing  the  steady  state  solution  is  as 
follows . 

Step  1.  Obtain  steady  state  solution  for  e[p]  using 
recursive  approximate  solution. 

Step  2.  Use  the  above  steady  state  solution  as  an  initial 
error  covariance  matrix  for  the  exact  solution 
(equation  (9))  and  propagate  for  several  time  steps. 

Step  3.  If  after  several  time  steps  the  exact  formulation 
shows  no  significant  change  in  expected  filter 
error  covariance,  the  steady  state  solution  of 
step  1 is  judged  valid. 

The  reasoning  for  the  above  test  is  that  if  the  steady  state  solution  of 

step  1 is  significantly  in  error,  a significant  change  in  expected  error 

covariance  should  show  up  in  step  3. 

Using  the  above  procedure  and  the  modeling  assumptions  discussed 

earlier  yields  the  results  shown  in  Table  B-1. 

Table  B-1.  Modeling  Assumptions 

Steady  State  Solution 
Propagated  7 Steps 

RAS  Steady  (7/16  sec)  using 

State  Solution  Exact  Formulation 

a = 1.68186  a = 1.68403 

r r 

= 0.43213  = 0.43248 

o =■  2.15636  o = 2.17497 

r r 

= 0.49593  CT^  = 0.49857 

I 


Probability  of 
Measurement 
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Notice  the  discrepancy,  indicated  by  last  column  of  Table  B-1,  is 
very  small.  As  in  the  transient  experiment,  expected  acceleration  error 
variance  is  virtually  identical  for  both  solutions  and  was  therefore  not 
tabulated. 

From  several  experiments,  such  as  the  two  presented  above,  it  was 
concluded  the  recursive  approximate  solution  is  quite  accurate  enough  for 
its  intended  purposes. 
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